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Abstract 

Raman microspectroscopy was used to characterize amorphous sihca plastic behavior. Using a 
correlation between Raman spectrum and density, a map of the local residual indentation-induced 
densification is obtained. The existence of a densification-induced hardening is also evidenced 
through a diamond anvil cell experiment. Such observations are not accounted for by the previously 
proposed hardening- free pressure-dependent yield criterion based on indentation curves. These 
results open the way to more accurate description of a constitutive law for amorphous silica. 
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I. INTRODUCTION 

Brittleness is a well-known and extensively studied characteristic of silicate glasses 
Comparatively, their plastic behavior has remained far less explored. Yet, the latter deter- 
mines the contact mechanics of the glass surface which is of primary interest for many 
micron-scale applications. 

For obvious historical reasons, the standard vision of plasticity is still largely focused on 
crystalline mater.als (especally metak). where irreversible stram .s tri^ered by shear stress 
and accommodated through volume-conserving dislocations motion [3J. Glasses do not fit 
into this frame for two different reasons. First, because of their amorphous nature there is 
no direct equivalent to these dislocation-based mechanisms Second, being amorphous, 
;heir structures, significantly more open than those of the chemically equivalent crystals 
^, makes them prone to densify. Indeed, silicate glasses densify ^\ while experiencing a 
plastic behavior that can be activated even by a purely hydrostatic loading This last 
characteristic leads to a distinction between normal and anomalous glasses, which exhibit 
respectively little and large densification, as well as different indentation fracture patterns 
0, Q- For fused silica (density 2.2 g.cm~^), the most anomalous silicate glass of industrial 
interest, this densification process is far from being negligible . It is acknowledged to 
densify up to about 20 % under mechanical stress |9|. Although this anomalous behavior of 
fused silica has been known for half a century in-depth investigations were limited by 
the brittle nature of silica, which restricts the characteristic size of plastic deformation to 
micrometer scale ^|. Indeed, evidences of the densification process at room-temperature 
can only be obtained experimentally in confined geometries, such as in a Diamond- Anvil Cell 
0] (DAC) or under an indenter jj, 3|, which do not allow direct access to local stress/strain 
data through the commonly used macroscopic techniques. 

We show here that, for amorphous silica, such local data of primary importance can 
be obtained experimentally through Raman microspectroscopy, which has recently proven 
particularly interesting as a method for local material characterization Q]. Through 
the combination of a specific method to obtain unaltered cross-sections of indentations with 
Raman microspectroscopic mappings, we provide a map of the local indent at ion- induced 
densification in a fused silica sample. This experimental data on the local strain state of the 
material constitutes a relevant sieve for the validation of a constitutive law for amorphous 
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The most sophisticated constitutive law for amorphous silica Q| does not account for 
these results. We suggest it be refined further by taking into account amorphous silica 
densificat ion-induced hardening, which we also evidence here with a DAC experiment. 

II. EXPERIMENTAL PROCEDURE 

A. Using Raman spectrum as a density gauge 

Between 200 and 750 cm~^, the Raman spectrum of silica glass consists in three bands 
(Fig. ^ contrasts the Raman spectrum of non-densified fused silica with the spectrum of 
largely indentation-densified silica glass). The main band at 440 cm^^ results from the 
symmetric stretching mode of bridging oxygens between two Si atoms. This band is intense, 
and affected by the densification process through a decrease of the inter-tetrahedral angles 
Si-O-Si 0], but it was shown to be also very sensitive to elastic loading j^l (such as 
the residual elastic strain expected in indentation). 

The defect lines Di and D2, at 492 and 605 cm~^, are respectively attributed to the 
breathing modes of the four-membered and three-membered rings |2^ . Their intensity ratio 



y,Q, 



was previously used in literature as an indicator for "fictive temperature" 
slight change of which corresponds to a slight change in density. Moreover, they shift in 
position towards higher wavenumbers with increasing density, as the associated breathing 
modes get frustrated by the lack of free volume. 

The D2 line is particularly useful because it has almost no overlap with the main band. 
Sugiura et al. correlated the position of the D2 line with the ratio of the sample density 
p to its initial density po- From their work, we get the following empirical relation: 

0.143 logio f-1 = \og,J-) (1) 



V and vq being respectively the final and the initial band position. The residual density 
evaluated through this relation accounts for both irreversible and elastic densification due 



to residual elastic strain. However, several works 



27l | evidenced that the D2 line position 



is only marginally sensitive to residual elastic strains. 
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B. Diamond- Anvil Cell geometry 



An amorphous silica (Saint-Gobain Quartz IDD) splinter with a ten-micron characteristic 
length was used as experimental sample. It was introduced into a DAG apparatus along with 
a very small piece of ruby. A Raman spectrum of the sample was acquired prior to the filling- 
up of the cell compartment with the pressure transmitting fluid. We used a transmitting 
iuid close to 1:5 methanol-ethanol mix, which insures hydrostaticity up to 16 to 20 GPa 
3|. That choice was motivated by the fact that the densification process in amorphous 
silica is expected to saturate between 20 and 25 GPa under hydrostatic loading j29i |. 

The sample was then loaded to a maximal applied hydrostatic pressure evaluated to 14 
GPa, using the shift of ruby fluorescence bands. After complete unloading, the Raman 
spectrum of the sample was measured. 



C. Indentation geometry 

Samples were prepared according to the following procedure. Plates of amorphous silica 
(IDD, Saint-Gobain Quartz), of dimension 10x10x0.4 cm^, were indented with an instru- 
mented microindentation device allowing surface visualization. These plain 2-kg-Vickers 
indents allowed us to gain access to top views of the densified area. 

2-kg-Vickers indents were also performed on a previously subcritically-grown crack, as 
described on figure |21 This provided us, on the crack surface, with a cross-section view of 
the densified material area located under the indent. In order to minimize the mechanical 
perturbation due to the presence of the crack, we performed a symmetrical loading. The 
indenter was thus positioned so that one of the Vickers pyramid diagonals be superimposed 
with the crack. Indentation curves obtained for on-crack indentations were similar (figure 
ini shows that both indentation curves almost coincide) with that obtained for plain 2-kg- 
Vickers indentations, showing that the presence of a crack did not noticeably affect the 
indentation process. 

We used a Raman microspectroscopy device (RMIOOO, Renishaw) including a 30 mW 
514 nm Argon laser, under a xlOO objective. The laser beam was focused on the surface 
without confocal system. We however checked on one section view that the mapping obtained 
using confocal system agreed with that obtained without confocal system. We evaluated the 
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excited volume of amorphous silica to approximately 2x2x5 /im^. The characteristic size 
of that scanned volume, which directly depends on the laser wavelength, is too large to 
map an elastoplastic nanoindent. However, it can easily be used on micro indents, such as 
that obtained with 2-kg-Vickers indents. Under the assumption that microcracks have no 
significant effect on the indentation-induced strain field, such experimental data can be used 
to validate an elastoplastic constitutive law for amorphous silica. 

After checking densification symmetry with respect to the indenter geometry, we per- 
formed surface mappings. Owing to the symmetry of the Vickers indenter, we first mapped 
one eighth of the densified area on three different plain indents. We thus recorded spectra 
in 13 locations with a 10 /im-spacing, according to Fig. HJ during 1500 s each. 

We then mapped half of the densified area on three different section views, recording 
spectra in 42 locations with a 5 yum-spacing, according to Fig. We also performed 
similar measurements on subsurface positions in order to complete the section map. These 
measurements are to be compared with those obtained under indent edges on top view 
mappings. For each spectrum, Raman scattering signal was accumulated during 900 s. 

III. RESULTS 

A. Diamond-anvil cell geometry 

Figure ini contrasts the initial spectrum of the sample (plain line) and that obtained after 
a 14 GPa hydrostatic loading (dotted bold line). Using relation (^J) between the D2 line 
shift and densification, the densification experienced by the sample is derived directly from 
the spectra measured before and after hydrostatic loading. Thus the sample experienced a 
densification of approximately 10%. 

B. Indentation geometry 

The D2 peak shift was then used as a density estimator according to equation (Q). A 
typical top view map is displayed on Fig. HI A typical cross-section map is displayed on 
Fig. El 

The iso-densification lines obtained for the top view map (Fig. are concentric. Their 
shape does account for the anisotropy of the Vickers pyramid. Indeed, densification is 
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more pronounced under the indenter edge than it is under the indenter face. This result is 



consistent with the fact that shear stress 
Vickers pyramid) enhances densification 



necessarily more intense under acuter parts of the 



30|. 



The iso-densification lines obtained for the cross-section map (Fig. ^ are shaped into 
"concentric bowls" , which agrees with the optical observations by Hagan and Van der Zwaag 
0]. Quantitatively, we evaluated the densification just under the indenter tip at around 
16 %. These results are consistent with the saturation of the densification process that is 
accepted to occur around 20% . Moreover, we notice that the subsurface measurements are 
consistent with those obtained along the indent diagonal on top view maps. This validates 
that the indents performed on pre-existing cracks are indeed reasonably equivalent to plain 
indents. The densification just under the indenter tip was again estimated to approximately 
16 %, in agreement with the previous measurement. 

IV. DISCUSSION 

A. Yield criterion for amorphous silica 

In most crystalline materials, plasticity occurs when the shear component of the applied 
stress reaches a yield value. Formally, this yield criterion can be written : 

fig) = r~Y (2) 

where r represents the intensity of shear stress and Y the yield strength. As long as 
fin.) remains negative (i.e. as long as the shear stress remains smaller than the material 
yield strength), the material behaves elastically, while it may behave plastically only when 
f{a) = 0. However, such a criterion cannot account for the silica glass densification under 
purely hydrostatic loading Q]. In order to do so, a dependance of /(£) to the hydrostatic 
pressure p has to be taken into account. 

Recently, a hydrostatic-pressure-dependent yield criterion was proposed to describe the 
plastic behavior of dense amorphous materials like silica glass 
yield criterion, /(a), can be written: 



was proposed to 

0,Q|iiQ,Q. 



This plastic 



(y S ' S 

fig) = ap+{l-a)r-Y = --Tr{g) + (1 - «)^^ - Y (3) 
where Tr{g) = an + + o"33 and g = g — Tr{g^l_ is the deviatoric part of 
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This pressure-dependent yield criterion is parameterized by two physical characteristics. 
One of them is the yield strength, Y. The other one, known as the densification factor 
[3] balances the respective influence of applied hydrostatic pressure, p, and shear stress, 
r.For a = 0, we obtain the Van Mises criterion (Eq. ((21) ), which describes the densificat ion- 
free plasticity of metals and, for a = 1, we describe a densification-only plastic behavior 
triggered exclusively by hydrostatic pressure, provided the plastic behavior is associated. 
Hence the denomination of a as the densification factor: the greater a, the more anomalous 
the material. 

Under the assumption of an associated perfect-plastic behavior, Xin and Lambropoulos 
[isj l adjusted these two parameters by comparing experimental load/displacement indenta- 
tion curve with finite element simulations. With the values a = 0.6 and Y = 5.43 GPa, 
they reproduced their experimental data: a = 0.6 accounts for a large contribution of den- 
sification to plastic strain (even allowing fused silica to exhibit a slight contraction of the 
cross section while being uniaxially compressed); and the hydrostatic loading yield pressure, 
obtained for r = 0, is consistent with DAC experimental results {2^. 



B. Amorphous silica densification-induced hardening 

We provide here a simple argument proving that amorphous silica hardens. Let us assume 
amorphous silica has a perfect plastic behavior (that is no hardening), at least until it 
reaches saturation density (See figure [7j). When hydrostatically loaded, before saturation 
occurs, the applied stress on amorphous silica will never exceed the material yield strength, 
which is the only applied stress at which plastic flow occurs in a perfect plastic material. 
Thus, amorphous silica has necessarily reached saturation density before it can reach an 
applied pressure higher than its yield stress. Our DAC sample was loaded above 9 GPa, 
which is the yield pressure provided by Xin and Lambropoulos' model. However, its 
densification is evaluated to a mere 10%, which is lower than both literature claims (20%) 
and our measurements under the indenter tip (16%) for the saturation density. This leads 
to the conclusion that amorphous silica indeed exhibits hardening. Moreover, similar DAC 
cycling-loading experiments would provide us with the evolution of amorphous silica yield 
criterion with material density along the hydrostatic pressure axis. 
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C. Densification map: a relevant local characterization of indentation-induced 
strain field 

Indentation experiments result in largely inhomogeneous stress distributions. The applied 
force (resp. the penetration) is merely the integral of local applied pressure (resp. vertical 
displacement) and as such simply proportional to the mean applied pressure (resp. vertical 
displacement). The exclusive use of such integral data may prove insufficient to identify the 
parameters of a constitutive law trying to account for more local strain measurements. Our 
Raman microspectroscopic maps may constitute a severe validation for such laws and it may 
help to gain insight on the refinements to perform. 

To test the infiuence of the densification factor on the general shape of the residual 
densified area, we performed FE calculations, using the software Abaqus, to model a cone 
indentation (of semi-apical angle of 70.3° — making it the "axisymmetric equivalent" to a 
Vickers indenter). Consistently with Xin and Lambropoulos' approach the material 
plastic behavior was assumed perfect and associated. The saturation, accounted for by pre- 
venting the material from densifying or fiowing any longer, occurs when material permanent 
densification reaches 16%. Calculations were performed for several values of a ranging be- 
tween 0.2 and 0.6. The yield strength, F, was chosen to match the experimental value in 
hydrostatic conditions {Xjoi = 9.05 GPa). Residual densification maps for a = 0.2 and 
a = 0.6 are presented in Fig. |H| When comparing these two maps, we notice that the 
indentation-densified areas are markedly different. In the case a = 0.6, with a large depen- 
dence on hydrostatic contribution, the numerical results predict a small but almost entirely 
saturated densified area . However, in the less anomalous case a = 0.2, the densification area 
is larger but with a smoother transition towards the non-densified material. This contrast 
highlights the direct relevance of the density map for the identification of the constitutive 
law. 



D. The role of strain- hardening in amorphous silica plasticity 

Let us now compare our experimental density map with the numerical results proposed by 
the literature. The densified area obtained by Xin and Lambropoulos j3l with a = 0.6 and 
Y = 5.43 GPa is qualitatively close to that provided by figure |Hfb). The general shape and 
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extension of the densified area is recovered. This confirms that the densification contribution 
to the sihca plastic behavior is significant Il5l |. However, the mainly saturated core 

with a sharp near-boundary density gradient largely contrasts with the more progressive 
densification of the experimental result. This shows clearly that the constitutive law for 
silica densification identified using indentation curves needs to be further refined to account 
for the indent at ion- induced densification map. 

At least one of the assumptions made in the constitutive law has to be relaxed. We 
expect that taking into account the densification-hardening process, the existence of which 
we have demonstrated by DAC measurements, will result in the concentric "bowl-shaped" 
areas associated to the experimental smoother density gradient, without much change in the 
global dimensions of the densified area. 



V. CONCLUSION 



Raman microspectroscopy provides experimentally relevant local data on the indentation- 
induced densification in amorphous silica. They constitute a discriminant sieve for the 
identification of its constitutive law, as they cannot be accounted for by existing models 
identified through indentation curves. Our results also clearly evidence a densification- 
induced hardening process. Taking it into account would allow one to model more accurately 
the plastic behavior of amorphous silica. 

We are currently working on finite-element simulation to provide such a plastic criterion 
for amorphous silica behavior. On-going experimental work also investigates the use of 
microspectroscopic methods to provide useful information on the mechanical behavior of 
other more complex silicate glasses. 
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Figure Captions 

Figure 1: Raman spectra for amorphous silica. The plain line spectrum was obtained 
with 2.2 g.cm"'^ amorphous silica, while the dotted line corresponds to indentation-densified 
amorphous silica. 

Figure 2: Experimental procedure allowing the obtention of a section view of the in- 
dented area. Using a scratch as surface defect (a), a subcritical crack is grown in the sample, 
using a double-bending loading geometry (b). Symmetrical indentation is then performed 
on the crack (c). Crack propagation is then completed under double-bending (d). Two twin 
sections of the indented area are then obtained. 

Figure 3: Indentation curves obtained for both on-crack indentation and plain inden- 
tation. They almost superimpose which shows us that the presence of a crack did not 
significantly affect the indentation process. 

Figure 4: Densification map of one eighth of a top view of an indentation-densified area 
obtained using the D2 band position indicator. The center of the indent is at coordinates 
(0;0). The bottom thin line figures the edge of the indent while the right thin line cor- 
responds to the middle of the pyramid face. The crosses correspond to the measurement 
points. The bold lines are the iso-densification lines. Numerical values are given in percent. 
Iso-densification lines appear to be concentric and account for the indenter anisotropy, as 
densification appears more important under the indenter edge. 

Figure 5: Densification map of half a section view of an indentation-densified area ob- 
tained using the D2 band position indicator. The center of the indent is at coordinates (0;0). 
The top bold line represents the sample surface. The crosses correspond to the measure- 
ment points. The black lines are the iso-densification lines. Numerical values are given in 
percent. Iso-densification lines appear to be shaped into concentric bowls in agreement with 
observations by Hagan et al. 

Figure 6: Raman spectrum obtained with an amorphous silica sample before (plain line) 
and after (dotted line) a 14 GPa hydrostatic loading. 

Figure 7:Expected behavior of a perfect plastic material under hydrostatic loading. Be- 
fore pressure reaches the yield pressure {pq = 9 GPa) the sample does not densify. As its 
plastic behavior is hardening-free, applied pressure cannot reach higher values than Pq until 
irreversible densification, AV/V, reaches its saturation value AV/V\sat- 
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Figure 8: Residual densification maps obtained through FE simulation of a 2D- 
axisymmetric cone indentation. They were obtained for two different values of the den- 



sification factor, (a) a = 0.2 and ( 
results by Xin and Lambropoulos 



3) a = 0.6. Densification values are given in percent. The 
15l | are similar to (b) 
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FIG. 1: Raman spectra for amorphous silica. The plain line spectrum was obtained with 2.2 g.cm 
amorphous silica, while the dotted line corresponds to indentation-densified amorphous silica. 




FIG. 2: Experimental procedure allowing the obtention of a section view of the indented area. 
Using a scratch as surface defect (a), a subcritical crack is grown in the sample, using a double- 
bending loading geometry (b). Symmetrical indentation is then performed on the crack (c). Crack 
propagation is then completed under double-bending (d). Two twin sections of the indented area 
are then obtained. 
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FIG. 3: Indentation curves obtained for both on-crack indentation and plain indentation. They 
almost superimpose which shows us that the presence of a crack did not significantly affect the 
indentation process. 
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FIG. 4: Densification map of one eighth of a top view of an indentation-densified area obtained 
using the D2 band position indicator. The center of the indent is at coordinates (0;0). The bottom 
thin hne figures the edge of the indent while the right thin Une corresponds to the middle of the 
pyramid face. The crosses correspond to the measurement points. The bold lines are the iso- 
densification lines. Numerical values are given in percent. Iso-densification lines appear to be 
concentric and account for the indenter anisotropy, as densification appears more important under 
the indenter edge. 
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FIG. 5: Densification map of half a section view of an indentation-densified area obtained using 
the D2 band position indicator. The center of the indent is at coordinates (0;0). The top bold line 
represents the sample surface. The crosses correspond to the measurement points. The black lines 
are the iso-densification lines. Numerical values are given in percent. Iso-densification lines appear 
to be shaped into concentric bowls in agreement with observations by Hagan et al. 
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FIG. 6: Raman spectrum obtained with an amorphous silica sample before (plain line) and after 
(dotted line) a 14 GPa hydrostatic loading. 
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FIG. 7: Expected behavior of a perfect plastic material under hydrostatic loading. Before pressure 
reaches the yield pressure {p^ = 9GPa) the sample does not densify. As its plastic behavior is 
hardening-free, applied pressure cannot reach higher values than Pq until irreversible densification, 
AV/V, reaches its saturation value AV/V\sat- 
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FIG. 8: Residual densification maps obtained through FE simulation of a 2D-axisymmetric cone 
indentation. They were obtained for two different values of the densification factor, (a) a = 0.2 and 
(b) a = 0.6. Densification values are given in percent. The results by Xin and Lambropoulos [isl 
are similar to (b) 
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